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ABSTRACT 

While satellite galaxies within groups and clusters are more likely to exhibit quiescent 
star formation rates (SFR) than central ('field') galaxies, an enhanced quiescent like- 
lihood extends to galaxies that are several virial radii beyond groups/clusters. Using 
a galaxy group/cluster catalog from the Sloan Digital Sky Survey Data Release 7 to 
examine galaxies of M star > 10 9 ' 7 Mq, we show that this enhancement has two causes. 
First, massive halos are clustered, so neighboring halos that host their own satellite 
galaxies boost the observed quiescent fraction out to large distances. Second, and more 
interestingly, central galaxies exhibit a strongly enhanced quiescent fraction out to 2.5 
virial radii (with some enhancement out to 5 virial radii) beyond groups/clusters. We 
show that this enhancement for central galaxies can be explained by 'ejected' satel- 
lite galaxies that orbit beyond their host halo's virial radius. Using a cosmological 
iV-body simulation, we examine the orbital histories of ejected satellites, showing that 
they extend out to 2.5 virial radii beyond their host halo and compose up to 40% of 
all central galaxies out to this distance. After ejection, these satellites typically orbit 
for several Gyr beyond the virial radius before falling back in, during which time they 
(continue to) lose significant halo mass. We show that a model in which ejected satel- 
lites experience the same SFR evolution and quenching as satellites within a host halo 
can explain essentially all environmental dependence of galaxy quenching beyond a 
host halo's virial radius. The strong halo mass stripping of ejected satellites also re- 
duces the average halo masses of galaxies near massive groups/clusters, an effect that 
is potentially observable via weak lensing. Overall, the SFRs/colors and stellar/halo 
masses of ejected satellites present challenges for the standard halo model for galaxy 
occupation and highlight the importance of environmental history in galaxy evolution. 

Key words: methods:numerical - galaxies: clusters: general - galaxies: evolution - 
galaxies: groups: general - galaxies: haloes - galaxies: star formation. 
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1 INTRODUCTION 

Environment impacts galaxy evolution. Galaxies in groups 
and clusters are more likely to have suppressed ('quiescent') 
star formation rates (SFR) and lie on the red sequence, have 
more elliptical/bulge-dominated morphologies, and have less 
atomic/molecular gas than galaxies of similar stellar mass 
in comparatively isolated 'field' environ ments (|Oemler 
• Dressler|[l980l |Dressler fc Gunn|[l983l |Balogh et aL] 
Blanton fc Moustakas| ~ 



1974 



1997 



2009| for review). While these dif- 



ferences are strongest for galaxies in massive clusters, they 
persist in groups at least down to virial masses of ~ 10 12 M Q 
(e.g., Weinmann et al.|200 6 Wctze let al.|2012 l, comparable 
to the Milky Way and M31, whose satellites also show such 
environmental trends (Mateo 1998, for review). 

Several observational works have demonstrated that 



this environmental dependence is confined mostly to 
physical scales corresponding to the virial radius of a 
group/cluster and is driven by satellite galaxies, which are 
not the massive 'central' galaxy at the core of a host halo 
(e.g., |Hogg et al. 2004 |Kauffmann et al 



2008; Tinker et al. 



2011 



2004 



Weinmann 



et al.||2006| |Blanton fc Berrind||2007[ |van den Bosch et al 



Wetzel et al. 



20121. These re- 



sults are physically motivated, given that a host halo's virial 
radius broadly corresponds to a physical transition from 
the low-density 'field' environment to a high-density region 
where dark matter and gas are virialized. 

While galaxies within groups/clusters strongly exhibit 
the above environmental dependencies, such trends can ex- 
tend to galaxies beyond the virial radius of a group/cluster. 
Several authors have noted an enhanced quiescent/red frac- 
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tion and reduced HI emission in galaxies out to 2 — 4 R v 



2000 



Solanes et al. 



around massive clusters (Balogh ct al 
20(32lJHansen et al.|2009||von der Linden et ai".|2010||Wetzel 
et al.|2012| ), as well as for dwarf galaxies around lower mass 



groups (Y. Wang et al. 2009 Geha et al. 20121, including 



those near the Milky Way and M31 ( jFraternali et al.|2009 1 



The above observational trends are plausibly driven, 
at least in part, by these galaxies having passed within 



much smaller distances from a group/cluster. Indeed, Balogh 



et al. (20001 first noted from iV-body simulations of clus- 
ters that particles that have passed within the virial radius 
can then orbit back well outside it. Satellite galaxies can 
do so as well, as investigated both through analytic argu- 



ments (Mamon et al. 20041 and in numerical simulations 



(|Gill et al.|2005 
low et al.||20Q9| 


|Sales et al.|2007|| 
Wang et al. 112009 


Warnick ct al. 2008; |Lud- 
"[Knebe et al-pOllJ |Ma- 


hajan et al. 


2011 


Teyssier et al. 


2012| Bahe et al.| 2013 1. 



In general, these works found that these 'ejected' (or 'back- 
splash') satellites typically orbit back out to a maximum 
distance of 4i?200c ~ 2.5 -R20011Q beyond a host halo after 
passing through it, and that as many as half of all galaxies 
within this distance are composed of these ejected satellites, 
with higher fractions for less massive galaxies and around 
more massive host halos. Thus, ejected satellites are po- 
tentially critical in understanding the properties of galaxies 
near groups/clusters and in obtaining a complete picture of 
environmental dependence. 

In this paper, we explore observationally and theo- 
retically the evolution of star formation in galaxies near 
(1 — 7J?200m) groups and clusters, to understand in detail the 
physical extent of the environmental dependence of galaxy 
evolution and the importance of ejected satellites. We use 
a galaxy group/cluster catalog constructed from the Sloan 
Digital Sky Survey (SDSS) together with a cosmological N- 
body simulation, to which we apply the same group-finding 
algorithm for robust comparison of model results to observa- 
tions. This combination allows us to pursue a more rigorous 
and statistically significant investigation, including account- 
ing for the importance of correlated structure (neighboring 
halos). We test the importance of star formation quench- 
ing in ejected satellites by employing the same empirically 
motivated and calibrated model for satellite SFR evolution 



that we developed in Wetzel et al. J2012kb. To summarize 



our main result, we find that the quenching of star forma- 
tion in ejected satellites is consistent with proceeding in the 
same manner as for satellites that remain within their host 
halo, and therefore, ejected satellites can account for essen- 
tially all of the observed environmental dependence of galaxy 
quenching beyond the virial radius of massive host halos. 

This paper represents the fourth in a series of five. 
In |Tinker et al.| ( |2011[ ), hereafter |Paper i] we described 
our SDSS galaxy sample, presented our method for obser- 
vationally identifying galaxy groups/clusters, and showed 
that galaxy quiescent fractions are essentially independent 
of large-scale environment beyond their host halo. In | Wetzel] 
et al. (2012 |a, hereafter Paper II we used our SDSS group 



1 We define a halo's virial radius, such that the aver- 

age interior density is A times a reference density, p re f: = 
A±irp Ie[ R 3 A . A = 200m means 200 X the average matter density 
of the Universe, and A = 200c means 200 X the critical density. 



catalog to examine in detail the SFR distribution of satel- 
lite galaxies and its dependence on stellar mass, host halo 
mass, and halo-centric distance, finding that the SFR dis- 
tribution is bimodal in all regimes. In Wetzel et al. (2012 lb, 



hereafter jPaper lJI\ we combined the above results with a 
cosmological iV-body simulation to constrain the star forma- 
tion histories and quenching timescales of satellites, showing 
that their star formation must evolve in the same manner as 
central galaxies for 2 — 4 Gyr (depending on stellar mass) af- 
ter falling into a host halo, but that once quenching starts, it 
occurs rapidly, with a SFR fading e-folding time of < 1 Gyr. 



Finally, in Wetzel et al. (20131, hereafter Paper V we will 



use the detailed orbital histories from our simulation to test 
physical mechanisms responsible for satellite quenching. 

Before proceeding, we define some nomenclature. We 
refer to galaxies with low SFR (see [2.1.1 1 as being 'quies- 
cent', and 'quenching' is the physical process of SFR fading 
rapidly from actively star- forming to quiescent, under the 
ansatz that once a satellite is quenched it remains so indefi- 
nitely. For our galaxy group catalog, 'group' means any set 
of galaxies that occupy a single host halo, regardless of its 
virial mass, and we use 'host halo' as a more general term 
for group or cluster. Finally, we cite all masses using h — 0.7 
for the Hubble parameter. 



2 METHODS 

Here, we briefly describe our observational galaxy sample 
and cosmological simulation. For full details, see Papers I 
and II for the former and |Paper III| for the latter. 



2.1 



2.1.1 



Observational catalogs from SDSS 

Galaxy catalog from SDSS 



Our galaxy sample is based on the NYU Value-Added 
Galaxy Catalog (Blanton et al. 2005[ ) from Data Release 
7 of SDSS ( |Abazajian et al.||2009|| ."Stell ar masses are from 
the kcorrect code of |Blanton fc Roweis| ( |2007) > , assuming a 
Chabrier (2003) initial mass function (IMF). We construct 
two volume-limited samples of galaxies with M r — 5 \og(h) < 
— 18 and —19, which go out to z = 0.04 and 0.06, from 
which we identify stellar mass complete limits of 5 x 10 9 and 
1.3 x 10 10 Mq, respectively. We use specific star formation 
rates, SSFR = SFR/M star , based on the current releas^jof 
the spectral reductions of Brinchmann et al. (20041. We de- 



fine galaxies with SSFR < 10 Mq as 'quiescent', based on 
the bimodal nature of the SSFR distribution (see |Paper II[ ). 



2.1.2 Galaxy group catalog from SDSS 

We identify groups of galaxies that occupy the same host 
halo using a modified implementation of the group-finding 
algorithm of |Yang et al.| ([2005 2007). For our group catalog, 
we define host halos as spherical overdensities containing 
200 times the average matter density of the Universe (A = 
200m). As Paper I outlined, we assign halo virial masses to 
groups by matching the abundance of halos above a given 



http : //www .mpa-garching.mpg . de/SDSS/DR7/ 
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dark matter mass to the abundance of groups above a given 
total stellar mass, n(> M v ir,haio) = n(> M atar , group), using 



the host halo mass function from Tinker et al. ( 2008 1 based 
on a flat, ACDM cosmology of fi m = 0.27, ft b = 0.045, h = 
0.7, n s = 0.95 and as — 0.82. The center of a group is given 
by its most massive galaxy, which we call the 'central' galaxy. 
Every group contains one central galaxy and can contain any 
number (including zero) of less massive 'satellite' galaxies. 



measure the orbital distance and velocity for each satellite 
in a continuous manner by numerically integrating its orbit 
between outputs, assuming energy and angular momentum 
conservation in a spherical NFW ( Navarro et al.|[l997 l po- 
tential given by the host halo's mass and concentration. 



2.2 Simulation 

2. 2. 1 Simulation properties 

To track the evolution of satellites from first infall to fi- 
nal merging/disruption across a broad range of host halo 
masses, we employ a dissipationless, iV-body simulation, us- 



ing the TreePM code of |White| |2002j with flat, ACDM cos- 
mology of Q m = 0.274, O b = 0.0457, h = 0.7, n = 0.95 and 
cr 8 = 0.8, evolving 2048 3 particles of mass 1.98 x 10 8 M 
with a Plummer equivalent smoothing of 2.5/i -1 kpc in a 
250ft _1 Mpc box, and storing 45 outputs spaced evenly in 
ln(a) from z — 10 to 0. This simulation first was presented 



in White et al. (20101 



2.2.2 Halo and subhalo tracking 

In the simulation, we identify 'host halos' using the Friends- 
of-Friends (FoF) algorithm ( |Davis et al.|1985| with a linking 
length of b = 0.168 times the average inter-particle spacing, 
which groups particles bounded by an isodensity contour of 
~ 100 x the average matter density. Within host halos, we 
identify 'subhalos' as overdensities in phase space through a 
6-dimensional FoF algorithm (F0F6D). For both host halos 
and subhalos, we keep all objects with at least 50 particles, 
and we define the center and velocity via the most bound 
particle. We assign to each (sub)halo a 'child' (sub)halo at 
the next simulation output, based on its 20 most bound 
particles. We define a 'central' subhalo as being the most 
massive subhalo in a newly formed host halo, which almost 
always corresponds to the subhalo at the minimum of the 
potential well. A subhalo retains its central demarcation un- 
til it falls into (is linked by the FoF algorithm to) a more 
massive host halo, becoming a 'satellite' subhalo. Thus, ev- 
ery sufficiently bound halo hosts one central subhalo and 
can host any number of satellite subhalos. 

We assign to each subhalo a maximum mass, A/ max , as 
given by the maximum host halo mass that it ever had as 
a central subhalo, motivated by the expected correlation of 
this quantity with galaxy stellar mass (see below). For a 
central subhalo, M max almost always corresponds to its cur- 
rent halo mass, the primary exception being ejected satel- 
lites (see below). For a satellite, M max almost always corre- 
sponds to its halo mass shortly before first infall. Also, as 



demonstrated in Wetzel & White (20101, simulations at our 



resolution scale can resolve and track massive satellite sub- 
halos past the point at which the galactic stellar component 
that they host would (start to) be stripped, merge with the 
central galaxy, or otherwise be disrupted, and we use their 
scheme by removing satellites with M DO und /M max < 0.007. 

Finally, to ensure that the simulation's output time 
spacing of ~ 500 Myr at z < 1 does not bias measurements 
of satellite orbital properties, such as pericentric distance, we 



2.2.3 Defining satellite infall and ejection 

We define 'first infall' of a subhalo as when it first becomes 
a satellite in a more massive host halo, if it remains a satel- 
lite for at least two consecutive outputs. The latter criterion 
avoids cases of temporary bridging if a subhalo briefly passes 
through the outskirts of a host halo. After experiencing first 
infall, if a satellite orbits outside of its host halo, becoming 
(again) a central subhalo in a distinct halo, we consider it to 
be an 'ejected satellite'. Thus, an ejected satellite appears 
to be a central subhalo in an instantaneous sense, though as 
we will demonstrate, it is more appropriate to continue to 
consider it a satellite. Finally, if an ejected satellite grows 
in halo mass by > 50% after ejection, we define it to be a 
newly formed halo and discard it from the ejected popula- 
tion, though as we will show, this criterion affects < 2% of 
ejected satellites, regardless of satellite mass. 



2.2.4 Galaxy group catalog from simulation 

Under the ansatz that a galaxy resides at the center of each 
dark matter subhalo, we assign stellar mass using subhalo 



abundance matching (SHAM; Vale & Ostrikcr 2006 Con- 
roy et aL||2006 1, a method that assumes a one-to-one map- 



ping that preserves rank ordering between subhalo M, 
and galaxy M star , such that n(> M max ) = n(> Mstar). 
Here, one assigns M star to subhalos empirically through an 
observed stellar mass function (SMF) that is recovered, by 
design. SHAM has succeeded in reproducing many observed 
galaxy statistics, including spatial clustering, satellite frac- 
tions, cluster luminosity functions, and luminosity- velocity 
relations ( e.g., |Conroy et al.|2006| |Berrier et al.|2006| |Yang 



et al.|2009| I Wetzel fc White|2010||Trujillo-Gomez et al.|2011 



Reddick et al. 



(2009) 



2012) 



We use the SMF from |Li fc White 
based on the same SDSS NYU-VAGC sample as our 
catalog, including the same if -correction and IMF. We ap- 
ply SHAM to the simulation at z = 0.05, implementing a 



as suggested by many observations (e.g., Yang et al. 


2008 


More et al.||2009| |Wetzel & White||2010| Leauthaud et al. 


2012 


Reddick et al.|2012 


Cacciato et al.|2013 


1. 



To make robust comparisons with the SDSS group cat- 
alog, we produce a 'simulation mock group catalog' by ap- 
plying the same group-finding algorithm to the simulation. 
While we base our models of SFR evolution on true satel- 
lite/central demarcation in the simulation, we effectively 
'observe' our model results through this group catalog, us- 
ing the same host halo definition and including observational 
effects, such as redshift-space distortions, that affect group 
purity/completeness as well as satellite/central assignment. 
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Figure 1. From SDSS: fraction of galaxies that are quiescent 
(SSFR < 10 -11 yr^" 1 ) versus projected distance to center of 
cluster-mass host halos, for galaxies that lie within line-of-sight 
Av = ±2 V200m of the host halos. Panels show bins of galaxy 
stellar mass, while different curves show galaxies decomposed by 
type: all (black), satellites (both within the cluster-mass host 
halos and within neighboring host halos; red), and neighboring 
central galaxies (blue). Solid curves show average value in the 
distance bin, while dotted lines show average across the entire 
SDSS sample. Error bars show 68% confidence interval of the 
average for a beta distribution. Considering all galaxies, an en- 
hanced quiescent fraction extends out to 6i?200m) while central 
galaxies exhibit a strong enhancement out to 2.5-R200m- 



3 STAR FORMATION IN GALAXIES NEAR 
GROUPS AND CLUSTERS 

We first examine observationally the SFRs of galaxies in and 
near cluster-mass host halos. We improve our analysis from 
|Paper lT| in which we examined all galaxies regardless of 
their redshift at a projected distance from a host halo, by 
examining only galaxies that lie in close proximity to a given 
host halo via applying a line-of-sight velocity cut of Aw = 
±2 Vboom, in which V^oom = y^GMvoomJlhoom is the virial 
velocity. (For reference, V2oom(10 14 M ) = 550 kms -1 .) 

Fig. [I] shows the fraction of galaxies that are quiescent 
as a function of projected halo-centric distance, d pro j, in two 
different bins of stellar mass in/around cluster-mass host ha- 
los. Curves show the average quiescent fraction, stacking all 
host halos in the given range of M2oom, and error bars show 



68% confidence interval on the average, as given by a beta 
distribution (Cameron 20111. Considering all galaxies, re- 
gardless of type (solid black curves), those at d pro j < Rsoom 
have the highest likelihood of being quiescent, but a clear 
enhancement, beyond the cosmic average of all galaxies in 
the stellar mass bin (dotted lines), persists out to ~ 6 7?200m 
(see also Hansen et al.|2009 von der Linden eTaLpOlOl ) . 

What causes this enhancement out to such large dis- 
tances, well beyond the turn-around (maximum) distance 
from which matter accretes onto these clusters? To address 
this, we use our group catalog to decompose these galaxies 
physically by type: satellite (solid red curves) and central 
(solid blue curves). (Note that some central galaxies appear 
within i?2oom as viewed in projection.) This decomposition 
reveals two trends that cause the enhanced quiescent frac- 
tion for all galaxies out to 6i?200m. First, structure is cor- 
related: massive host halos are more likely to be found near 
each otherJ^Thus, the likelihood that a galaxy is in another 
massive host halo decreases with distance from a cluster, so 
the fraction of galaxies that are satellites, and the average 
host halo mass of those that are, decreases with distance. 
Because the quiescent fraction is higher for satellites than 
central galaxies, and for satellites it is higher in more massive 



host halos ( Paper II I , correlated structure can cause a gradi- 



ent in the quiescent fraction with distance even if there were 
no gradients for central and satellite galaxies separately. 

However, Fig. [I] also shows another, more interest- 
ing trend: central galaxies, which reside in their own host 
halo, exhibit a strongly enhanced quiescent likelihood out 
to ~ 2.5i?200m- While this enhancement is strongest around 
cluster-mass host halos, central galaxies near lower mass 
host halos show similar behavior, and a non-zero enhance- 
ment beyond the cosmic average for central galaxies extends 
out to w 57?200m in essentially all cases (see Fig. [7] in jSJ. 

What causes this enhanced quiescent fraction in central 
galaxies near massive host halos? It potentially could be an 
artifact of interloping galaxies caused by redshift-space dis- 
tortions or some other aspect of the group-finding algorithm. 
However, as we will show using our simulation mock group 
catalog in ij5] such observational effects cannot account for 
this strong excess (in agreement with Y. Wang et al.|2 009). 
This leaves two feasible physical causes: 



(a) Strong environmental processes might extend well be- 
yond the (spherical) virial radius of host halos, quenching 
central galaxies prior to their first crossing inside. 

(b) The observed enhancement could originate from 'or- 
dinary' environmental quenching within a host halo's virial 
radius and propagate to larger distances via 'ejected' satel- 
lites, which orbit beyond the virial radius after falling inside, 
becoming (once again) central galaxies in a distinct halo. 

Given its simplicity, we postulate that (b) most natu- 
rally explains the observed trends, as we will subsequently 



examine and test. We will discuss (a) in [ 7.2.2 



3 For the cluster-mass host halos in Fig. ^ the median projected 
distance (given our velocity cut) to the center of the nearest host 

M© is d 



halo with M 20 0m > 10 14 > 13 > 12 
respectively. In real space 
median rf/i?200ro values are roughly twice as large 



proj/R200m = 4.7, 1.4, 0.6, 

as measured in our simulation, the 
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d / R 200m d / R 200m d / R 200m 

Figure 2. From simulation: fraction of galaxies that have been ejected versus (real-space) halo-centric distance at z = 0. Panels 
show bins of virial mass of the host halos on which these profiles are stacked, while different curves show bins of galaxy stellar mass 
(amplitudes decrease monotonically with stellar mass). Values at d > i?200m show the fraction of central galaxies that are ejected 
satellites, n e j ecte( j(<i)/ncen(c(), while values at d < R200m show the fraction of satellites that experienced an ejected phase and have since 
fallen back in, n c j cctc ^(d) /n sa t(<i). Solid curves include just satellites that were ejected from the host halos on which these profiles are 
centered, while dashed curves include all ejected satellites, including those ejected from neighboring host halos. Shaded regions show 
68% confidence interval of the average. Ejection is more common for less massive galaxies and in/around more massive host halos, being 
most significant at f — 2.5-f?200m, though extending to much larger distances because of satellites ejected from neighboring host halos. 



ORBITAL EVOLUTION OF EJECTED 
SATELLITES 



In |Paper IHl we showed that a significant fraction of satel- 
lites inside a host halo spent time in an ejected phase be- 
yond the virial radius before falling back in. In this section, 
we investigate in more detail the orbital histories of ejected 
satellites, building on previous work but focusing on orbital 
properties that are (potentially) relevant to quenching star 
formation. While we examine galaxies and host halos in the 
same mass ranges as in our SDSS catalog, in this section we 
use the halo catalog directly from our simulation (not from 
the mock group catalog) and measure halo-centric distance 
in real-space, d, in order to highlight physical trends without 
projection effects or redshift-space distortions. 



4.1 Fraction of galaxies that have been ejected 

We first examine the importance and radial extent of ejected 
satellites around massive host halos. Fig.[2]shows what frac- 
tion of galaxies at a given (real-space) distance around host 
halos are ejected. We show the average value from stack- 
ing galaxies around host halos of the given virial mass, and 
shaded regions shows the 68% confidence interval on this 
average. Values at d > i?200m show what fraction of central 
galaxies are ejected satellites, n e j ecte d(d)/n cen (d). For refer- 
ence and continuity, Fig. [2] also shows values at d < i?200m, 
which indicate what fraction of satellites experienced an 
ejected phase and fell back in, n c j octc d(d)/n sat (d). 

The solid curves show satellites that were ejected specif- 
ically from the host halos on which these profiles are cen- 
tered. Across all galaxy and host halo masses, a non-trivial 
population of ejected satellites extends to ~ 2.5 -R200m be- 
yond their host halo, in good agreement with previous sim- 
ulation results (|Gill et al.|2005|[Warnick et al.|2008||Ludlow 



likely caused by an orbital energy boost during a multi-body 
encounter within the host halo ( Sales et al.|2 007). However, 
this phenomenon primarily affects satellites with particu- 



et al. 2009 Wang et al.|2009 I. This maximum distance is ap- 
proximately as expected based on the turn-around distance 
of infalling matter, if it conserves energy as it orbits back 
to apocenter after infall (Ma mon et al.|2004 l. A small frac- 
tion of ejected satellites extend to somewhat larger distance, 



larly low mass that can be boosted more easily (Ludlow 
et al. |2009" |: the vast majority of satellites in our mass range 
have orbital apocenters that are similar to or smaller than 
their turn-around distance, implying that any such orbital 
energy boosts are mild for satellites that we examine. 

By contrast, the dashed curves in Fig. [2] show the 
ejected fraction regardless of which host halo the satellites 
passed through. Thus, it also includes satellites that were 
ejected from neighboring host halos. In examining this total 
fraction, a strong excess still persists out to 2.5 -ffeoom, but it 
declines only gradually before reaching the cosmic average 
(given in Fig.|3| at ~ 10 ifeoom- Thus, while ejected satellites 
are most common near massive host halos, they are cosmi- 
cally ubiquitous. Because separating ejected satellites based 
on which host halo they passed through observationally is 
highly non-trivial, if not impossible, one must consider this 
total population when comparing to observations. 

Fig. [3] shows more directly how the fraction of central 
galaxies that are ejected satellites varies with galaxy stel- 
lar mass. First, the solid curve shows what fraction of all 
central galaxies are ejected satellites. This fraction increases 
at lower stellar mass, tho ugh it is always < 10% down to 
10 9 M (in agreement with Wang et al.|2009 1. Thus, ejected 
satellites are of modest importance to the overall population 
of central galaxies. However, as Fig. [2]showed, ejected satel- 
lites are particularly common near massive host halos, and 
the dashed curves in Fig. [3] show the ejected fraction for 
central galaxies that are within 1 — 2.5 -R200m of host ha- 
los of the given virial mass. This ejected fraction strongly 
increases with host halo mass and decreases with galaxy 
mass. Most likely, these mass dependencies arise from the 
increased efficiency of dynamical friction 'braking' of or- 
bital velocity when the ratio of satellite to host halo mass is 
smaller ( |Boylan-Kolchin et al.|2008|| Jiang et al.|2008"| . Note 
that 40% of all central galaxies with M s t ar < 10 lu Mq near 
cluster-mass host halos are ejected satellites. 
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Figure 3. From simulation: fraction of central galaxies that 
are ejected satellites versus stellar mass at z = 0. Solid curve 
shows fraction with respect to all central galaxies, fi e j ecte d/"cen, 
while dashed curves show fraction within 2.5i?200m of host halos 
of the given virial mass, n c j cctcd (< 2.5 i?200m, M200m)/n C cn(< 
2.5 -f?200m) A'ftoOm)- Satellite ejection is more common for less 
massive galaxies and near more massive host halos. 
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Figure 4. From simulation: cumulative distribution of minimum 
halo-centric distance to which satellites at z = have orbited, in 
bins of host halo virial mass. Solid curves show satellites that 
currently are ejected (beyond -R200m)> while for reference, dashed 
curves show all satellites that currently are within i?200m- The 
median pericenter for ejected satellites is 0.25 — 0.5-R200m, be- 
ing monotonically smaller in more massive host halos. Ejected 
satellites almost always experience a single pericentric passage. 



4.2 Minimum pericentric distance 

We next examine how far into their host halo these ejected 
satellites orbited, in order to understand the feasibility of 
their star formation being affected. We measure a satel- 
lite's d/R,2oom at each simulation output and identify the 
minimum ratio experienced since first infall, using orbital 
integration as given in §2.2.2| Fig. [4] shows the cumulative 
distribution of minimum d/ifeoom to which different satellite 
populations at z = have orbited. Solid curves show ejected 



satellites, which currently are beyond i?200m, while for ref- 
erence, dashed curves show all satellites that currently are 
within Rwom- Across the entire sample, the median pericen- 
tric passage of ejected satellites is ~ l/3i?200m, sufficiently 
deep into the halo's potential well to make it feasible that the 
majority have experienced quenching process(es). As Fig. [4] 
shows, the median for all satellites within the virial radius 
is ~ 0.15i?200m, so by comparison, ejected satellites have 
not orbited especially deeply. Though, it may be unfair to 
compare ejected satellites to all satellites, given that almost 
half of all satellites currently lie within 1/3 ifeoom- If instead 
we compare the pericenter distribution of ejected satellites 
to those currently at 0.9—1 -ffeoom, they are almost identical 
(not shown). Thus, the majority of ejected satellites are not 
on orbits with particularly unusual pericenters, they simply 
have sufficient orbital energy to bring them back beyond the 
virial radius. 

We find no dependence of any of these distance distri- 
butions on satellite mass, but Fig.[4]shows satellites split by 
the virial mass of the host halo passed through. While the 
distributions are similar at d/i?200m < 0.1, satellites orbited 
to somewhat smaller d/J?200m in more massive host halos, 
with median values of varying from 0.25 — 0.5-R200m- This 
trend is likely driven, at least in part, by inf ailing satellites 
having smaller orbital pericenters and circularities in more 
massive host halos (jWctzcl 20TT]). 

Finally, the pericentric distances experienced by ejected 
satellites in Fig. [4] are almost always (~ 90%) the result of 
a single pericentric passage between first infall and ejec- 
tion, also noted previously ( |Gill et al.||2005[ ). The ~ 10% of 
ejected satellites that experienced multiple pericentric pas- 
sages typically are on highly radial orbits and experienced 
several passages in/out of the host halo. 



4.3 Orbital timescales 

Our primary goal is to test whether ejected satellites experi- 
ence the same quenching of star formation as 'normal' satel- 
lites that remain within a host halo. As |Paper III] showed, 
SFR in normal satellites evolves unaffected for 2 — 4 Gyr 
after first infall, depending on stellar mass, after which it is 
quenched rapidly, with an e-folding time of < 0.8 Gyr. Thus, 
we now examine the orbital timescales of ejected satellites 
in the context of these quenching times. 

Fig. [5] shows the distribution of time since first infall 
(dashed) and time since ejection (solid) for ejected satel- 
lites at z = 0. We find little dependence of these orbital 
times on the virial mass of the host halo that they passed 
through, as shown, and no dependence on galaxy mass (not 
shown). Across all ejected satellites, the median time since 
ejection is 1.4 Gyr and since first infall is 4.8 Gyr. The distri- 
bution of time since ejection falls off rapidly, meaning that 
relatively few satellites are on (nearly) unbound orbits that 
keep them beyond the virial radius for a long time. Instead, 
almost all ejected satellites remain bound to their host halo. 
For ejected satellites that are crossing back inside the virial 
radius at z — 0, we find that the median time spent in an 
ejected phase is 3.6 Gyr. 

We also find that the median amount of time that 
ejected satellites spent within their host halo is 2.9 Gyr 
across our sample, again with no significant dependence on 
galaxy or host halo mass. This time-as-satellite is similar 
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Figure 5. From simulation: distribution of time since first infall 
(dashed) and time since ejection (solid) for ejected satellites at 
2 = 0, in bins of the virial mass of the host halo that they passed 
through. The median time since first infall is 4.8 Gyr and since 
ejection is 1.4 Gyr. These orbital times do not depend significantly 
on either host halo mass or satellite mass and are comparable to 
satellite quenching timescales. 



to the virial crossing time (2 -ffeoom/ V^oom) at the typical 
redshift when these satellites fell in (z ~ 0.5), as expected 
given that the vast majority experien ced a single pericent ric 



passage, and in good agreement with Wang et al. (2009)]^] 



Thus, compared to the quenching delay timescales of 
satellites from Paper III , the typical time-as-satellite is sim- 
ilar, and the typical time since first infall is considerably 
longer. As such, most ejected satellites easily have the po- 
tential to be quenched in the same way as normal satellites. 

4.4 Evolution of halo mass 

Having shown that most ejected satellites orbited relatively 
deeply into their host halo and have been satellites for a suf- 
ficiently long time to be (at least potentially) quenched, we 
next explore the evolution of their (sub)halo mass after ejec- 
tion. If an ejected satellite's (sub)halo starts to accrete mass 
again after ejection, in the same way as non-ejected halos, 
this accretion feasibly could lead to renewed star formation, 
even if quenching already occurred. Alternately, if ejected 
satellites typically do not accrete mass, then they are likely 
to remain quenched in the same way as normal satellites. 

To demonstrate the evolution of halo mass after ejec- 
tion, Fig.|6]shows the cumulative distribution of the fraction 
of halo mass that an ejected satellite has at z = relative 
to what it had immediately after ejection. Different curves 
show satellites ejected from host halos of different virial 
mass. Across all host halos, only ~ 10% of ejected satellites 
gained halo mass since ejection. Instead, ejected satellites 
almost always continue to lose mass from tidal stripping 



4 While | Wang et al.| ( |2009[ l found a shorter timescale for more 
massive satellites, they conjectured that this arose from a higher 
liklihood of transient fly-by's for more massive satellites. Our 
more conservative selection criterion, which requires an ejected 
satellite to have remained within its host halo for at least two con- 
secutive simulation outputs, leads to no such mass dependence. 
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Figure 6. From simulation: cumulative distribution of the frac- 
tion of halo mass that an ejected satellite has at z = relative 
to what it had immediately after ejection, in bins of the virial 
mass of the host halo that it passed through. Ejected satellites 
experienced somewhat stronger stripping in/near more massive 
host halos, with median fractions ranging from 0.7 — 0.9. Only 
10 — 20% of ejected satellites have gained mass since ejection; 
most continue to experience halo mass stripping, with stronger 
stripping near more massive host halos. 



in the dynamically hot environment surrounding a massive 
host halo. Fig. [6] shows that stronger stripping occurs near 
more massive host halos, with the typical mass loss since 
ejection being 40% near clusters (in good agreement with 
~ Warnick et al.||2008l), but only 10% near 



Gill et aLj|2005 

M 20 om = lO 17 ^ 3 ' M . We find no significant dependence on 
satellite mass. This continued halo mass loss strongly sug- 
gest that the vast majority of ejected satellites do not restart 
star formation from renewed gas accretion after ejection, but 
rather continue to be affected and quenched after ejection 
in the same way as satellites within a host halo. These re- 
sults also have important implications for the M sta r — Afhaio 
relation of ejected satellites, as we will explore in Sj6] 

4.5 Summary 

To summarize this section, the orbits of ejected satellites ex- 
tend to « 2.5i?200m beyond their host halo, though because 
of neighboring host halos and correlated structure, an en- 
hanced total population extends out to ~ lOifeoom- Lower 
mass satellites and those associated with more massive host 
halos are more likely to experience an ejected phase. After 
first infall, these satellites almost always experience a single 
pericentric passage before ejection, typically to ~ 1/3 i?200m , 
which is similar to satellites that are just within i?200m- 
Thus, the typical time that ejected satellites spent within 
the host halo, ~ 3 Gyr, is approximately the virial crossing 
time, and almost all remain on bound orbits and re-enter 
the host halo after ~ 3.5 Gyr. For satellites that currently 
are ejected, their typical time since first becoming a satel- 
lite is ~ 5 Gyr, longer than the 2 — 4 Gyr quenching delay 
time of satellites, and most (90%) continue to lose halo mass 
after ejection. These results strongly suggest that ejected 
satellites quench and remain quenched in the same way as 
satellites inside a host halo, which we next test directly. 
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5 QUENCHING OF STAR FORMATION IN 
EJECTED SATELLITES 

In !|3j we postulated that ejected satellites could explain the 
observed enhancement in the quiescent fraction for central 
galaxies near massive host halos. We now test this scenario 
directly by examining two contrasting models for the evolu- 
tion of SFR in ejected satellites: 

(a) 'satellite-like' evolution, in which the SFRs of ejected 
satellites are quenched in the same way as normal satellites, 

(b) 'central-like' evolution, in which the SFRs of ejected 
satellites evolve in the same way as normal (non-ejected) 
central galaxies. 

To test 'satellite-like' evolution, we apply the same 
physically motivated, two-stage model for the evolution of 
SFR in satellites from |Paper III| To summarize, we start by 
obtaining statistically accurate initial SFRs for satellites at 
their time of first infall via an empirically based parametriza- 
tion for the evolution of central galaxy SFRs out to z = 1, 
combining our SDSS group catalog with data at higher z 



from the Cosmic Evolution Survey (COSMOS) ( Drory et al 



2009 1 and the All-wavelength Extended Groth strip Interna- 



tional Survey (AEGIS) ( Noeske et al. 2007 1 . If a satellite was 
actively star-forming at its time of first infall, tint, obtained 
via the simulation, we parametrize its subsequent evolution 
as follows. Its SFR fades gradually, in the same manner as 
central galaxies, across a 'delay' time, tq, delay After this, 
it starts to be quenched, and its SFR fades exponentially 
across an e-folding time, tq, f a de : 

SFR sat (t) = 

{SFR cen (t) t < t Q , start 

SFR C en(tQ,st a rt)eXp|-t^^l| t > tg, start 

(1) 



in which *q, start = iinf +tQ, delay, and 
SFR ccn (t) = SFR (t-t f )exp 



t-tf 



(2) 



with tf being the time of initial formation, which we take 
to be t(z — 3). We obtain SFRo and r ccn in narrow bins of 
stellar mass using the measured M s t ar and SFR of actively 
star forming central galaxies in the SDSS group catalog (as- 
suming 40% stellar mass loss through supernovae and stel- 



lar winds by z = 0). As Paper III showed, this empirically 
constrained model provides accurate SFR distributions and 
quiescent fractions for satellites at z — across both stellar 
and host halo mass, with tQ, delay and tq, fade that depend on 
satellite mass but not significantly on host halo mass. (As 
|Paper V] will show, this model also reproduces accurately 
the dependence of SFR versus halo-centric distance within 
host halos.) We use the same quenching timescales as |Pa-| 
|per III| t Qt delay = 3.4 Gyr and r Qt fad c = 0.8,0.6 Gyr for 
log (M s tar/ M ) = 9.7- 10.1, 10.1 - 10.5, respectively^ 



5 These timescales, and all results in this paper, are based on 
our 'fiducial' parametrization in |Paper III| for the initial quiescent 
fraction for satellites at 4; n f. Using our 'alternate' parametriza- 
tion, along with the associated quenching timescales, does not 
change any results in this paper significantly. 



For a contrasting null model, we consider 'central-like' 
evolution, in which ejected satellites have the same instan- 
taneous distribution of SFRs as central galaxies. Physically, 
this scenario implies either that SFR in ejected satellites is 
not affected environmentally, so it evolves in the same way 
as in normal central galaxies, or that SFR in ejected satel- 
lites is quenched environmentally, but it resumes to being 
the same as in normal central galaxies after ejection. To im- 
plement this model, we proceed as above, but we re-assign 
SFRs to ejected satellites empirically by drawing randomly 
from 'isolated' central galaxies of the same stellar mass in 
the SDSS group catalog that do not lie within 3 i?200m in 
projection of a more massive host halo, to minimize any 
possible contamination from ejected satellites. 

The two models above provide SFRs at z — for all 
satellites in the simulation. We assign SFRs to central galax- 
ies by drawing randomly from isolated central galaxies of the 
same stellar mass in the SDSS group catalog. Thus, we as- 
sign SFR to galaxies in the simulation according to their 
physical satellite versus central demarcation, but we then 
measure the results using the mock group catalog to com- 
pare robustly with SDSS. 

Fig.[7]shows the quiescent fraction for central-identified 
galaxies as a function of projected distance from host halos 
of the given mass in each panel, comparing the results of 
the two models above against the SDSS group catalog. As 
in Fig. [I] points show the average observed value in the 
distance bin from SDSS, and dotted lines show the average 
for all central galaxies of the given stellar mass in SDSS. 
Note that, as compared with Fig. [T] Fig. [7] extends to host 
halos of much lower mass, where the upturn near ^?200m is 
weaker, but the excess beyond the cosmic average remains 
strong out to 5i?200m. Each panel shows two bins of galaxy 
stellar mass; we find that galaxies of higher mass show less 
enhancement, as expected from Fig. [2] 

Curves show the two model results as measured through 
the simulation mock group catalog. Dashed curves show 
'central-type' evolution for ejected satellites, in which all 
central galaxies have the same SFR distribution, regardless 
of formation history. We note two important points. First, 
even when viewed through the mock group catalog, this null 
model recovers well the average quiescent fraction for cen- 
tral galaxies (dotted lines). Some upturn occurs at i?200m, 
caused by redshift-space distortions, but it is weak compared 
to the observed enhancement. Thus, we conclude that ob- 
servational effects related to the group finder alone cannot 
account for the strong enhancement in the quiescent frac- 
tion for central galaxies: it is a real physical effect. Second, 
because this model recovers the cosmic average, it underesti- 
mates the observed quiescent fraction even out to « 5 -fooom- 
Thus, the importance of quenching in ejected satellites is not 
confined just to within 2.5 -R200m of host halos; it is critical 
to understanding the entire population of central galaxies. 

By contrast, solid curves show 'satellite-type' evolution 
for ejected satellites, in which their SFR is quenched in the 
same way as normal satellites. This model provides good 
overall agreement for the upturn within 2.5i?200m- Just as 
importantly, it also correctly enhances the quiescent frac- 
tions at larger d pro j , because of satellites ejected from nearby 
host halos, as Fig. [2] showed. 

We emphasize that this model for 'satellite-type' SFR 
evolution is simply an extension of our empirically con- 
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Figure 7. Comparing models to SDSS: fraction of central galax- 
ies that are quiescent versus projected distance to massive host 
halos. Panels show bins of host halo mass, different colors show 
bins of galaxy stellar mass. From the SDSS group catalog, points 
show average value in the distance bin, and dotted lines show av- 
erage for all central galaxies in the catalog (as in Fig.JTJ. Curves 
show models as 'observed' through the simulation mock group 
catalog: if ejected satellites have the same SFR as isolated cen- 
tral galaxies (dashed), or if SFR in ejected satellites evolves and is 
quenched in the same way as other satellites, according to equa- 
tion Q (solid) . Overally, the latter model successfully can explain 
the enhanced quiescent fraction in central galaxies out to 5 i?200m ■ 
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Figure 8. From simulation: median ratio of a central galaxy's 
stellar mass, M atar , to its (instantaneous) halo virial mass, M^ & i , 
as a function of stellar mass at z = 0, for all central galax- 
ies (dashed) and just ejected satellites (solid), as obtained via 
abundance matching using M max - Because of strong stripping of 
(sub)halo mass after infall and ejection, ejected satellites have 
2.5 X higher ratio of M s tar to M^aio than for all central galaxies. 



strained model from |Paper III] — which provided excellent 
agreement for satellites that remain within host halos — to 
the ejected satellite population. Furthermore, our hypothesis 
that ejected satellites are quenched (and remain quenched) is 
supported physically by their orbital histories, as ij4]showed. 
Thus, given the good agreement of this model with SDSS in 
Fig. [2] we arrive at the following main conclusions: 

• SFR evolution and quenching in ejected satellites is con- 
sistent with proceeding in the same manner as in normal 
satellites that remain within their host halo. 

• This quenching of ejected satellites naturally can ex- 
plain the enhanced quiescent fraction for central galaxies 
out to 5 7?200m beyond host halos. 

• The success of this model leaves little-to-no room for en- 
vironmental processes that quench central galaxies beyond a 
host halo's virial radius, that is, essentially all environmen- 
tal quenching is consistent with originating within a host 
halo's virial radius. 



6 M s tar - M ha i RELATION FOR EJECTED 
SATELLITES 

In addition to star-formation quenching, in this last section 
we explore the implications of the mass evolution of ejected 
satellites. As §4.4| showed, because of the tidal stripping of 
(sub)halo mass that occurs after first infall and continues 
after ejection, ejected satellites have significantly reduced 
halo masses as compared with normal (non-ejected) central 
galaxies of the same stellar mass. However, because a galaxy 
is smaller and more tightly bound than its (sub)halo, ejected 
satellites likely do not experience significant stellar mass loss 
from tidal stripping during their (typically single) passage 
through the host halo (Knebe et al. 20111. Thus, central 



galaxies that are ejected satellites would have a higher ratio 
of stellar to halo mass than normal central galaxies. 
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Figure 9. From simulation: average halo mass for central galax- 
ies in the given stellar mass bin versus halo-centric distance at 
z = 0, normalized to the average for isolated central galaxies in 
the stellar mass bin. Solid curves show values in real-space from 
the simulation directly, while dashed curves show values in projec- 
tion measured using the simulation mock group catalog. Shaded 
regions show 68% confidence interval of the average. Top: Qui- 
escent central galaxies near cluster-mass host halos, where the 
signal is strongest. Bottom: All central galaxies near group- and 
cluster-mass host halos, where statistics are stronger. A signif- 
icant reduction in average halo mass, driven primarily by tidal 
stripping of ejected satellites, persists out to 2.5-R200m- 



To illustrate the significance of this effect for the entire 
population at z = 0, Fig. [8] shows the median M s tar/Afh a i 
ratio for all central galaxies in the simulation (dashed) and 
for just central galaxies that are ejected satellites (solid), as 
obtained through subhalo abundance matching using M max 



(i 2.2.4 1. Because of continued mass stripping, the halo mass 
of an ejected satellites is typically only ~ 40% that of the 
whole central galaxy population of the same stellar mass. In 
other words, ejected satellites systematically have a ~ 2.5 x 
higher ratio of M a t ar to Mhalo than all central galaxies, with 
a somewhat stronger fractional enhancement at lower mass. 

While Fig. [8] shows the offset in halo mass if one iso- 
lates the population of ejected satellites, in reality they 
are spatially intermixed with normal (first-infalling) cen- 
tral galaxies. Because of the intrinsic difficulty in separating 
these populations, observationally measuring this effect re- 



quires stacking and averaging both populations, leading to a 
weaker measured offset. The population of ejected satellites 
is most significant near the most massive host halos (Fig.|2|, 
where halo mass stripping also is strongest (Fig. j6j). Addi- 
tionally, ejected satellites are more likely to be quiescent 
than non-ejected central galaxies of the same stellar mass 
(Fig. 17b. Thus, quiescent, low-mass central galaxies near the 
most massive host halos should provide the strongest mea- 
surable reduction in average halo mass. 

To examine the potential observability of this effect, 
we use our 'satellite-like' evolution model from ij5] to as- 
sign SFRs to all galaxies in our simulation. We then stack 
all quiescent central galaxies near cluster-mass host halos 
and measure their average halo mass as a function of halo- 
centric distance. Fig. [9] (top) shows this average halo mass, 
normalized to the average halo mass of all isolated (do not 
lie within 3 ifeoom of a more massive host halo) central 
galaxies of the same stellar mass. We select galaxies with 
M s tar = 10 9 - 7 " 10 - 5 M to match our SDSS sample, but we 
note that this effect is even stronger at lower stellar mass 
because of a higher likelihood of being an ejected satellite. 

Solid blue curves show the true, underlying average as 
measured in real space, using the halo catalog directly from 
the simulation. At d « ifeoom, where ejected satellites are 
most common, the average halo mass of all central galaxies 
is 45% that of isolated central galaxies)^] This reduction re- 
mains particularly strong out to ~ 2.5i?200m, as expected 
given the spatial extent of ejected satellites from Fig. [2] 
Similar to Fig. [7J the average ratio does not asymptote un- 
til « 10 -ffeoom because of satellites ejected from neighbor- 
ing host halos, and because ejected satellites are cosmically 
ubiquitous, the ratio asymptotes to ~ 0.93 and not 1. 

Dashed red curves show the average as measured in pro- 
jection using the simulation mock group catalog, imposing 
a line-of-sight velocity cut of ±300 km/s (more conservative 
than before to maximize signal). While measuring in projec- 
tion reduces the signal somewhat, it remains strong if using 
a sufficiently conservative velocity cutQ 

These results suggest that the reduced halo masses for 
central galaxies near groups/clusters as caused by ejected 
satellites could be measured observationally via galaxy- 
galaxy weak gravitational lensing, a technique that can mea- 
sure accurate halo masses for samples of stacked galaxies 
(e.g.,[Mandclb aum et al.|20 06). However, in applying such a 
technique to the sample in Fig. [9] (top), note that stack- 
ing all clusters with M 20 0m = l(r 4 ~ 15 M Q in our SDSS 
group catalog leads to ~ 500 quiescent galaxies at M star = 



10 M 



Mq in each bin of Ad plo j/R2 



1, a number 



that is unlikely to yield a discriminating stacked lensing sig- 
nal from current or near-term surveys \Li et al.|2013[ ). How- 
ever, if we consider instead all central galaxies near host ha- 



6 Some tidal stripping does occur as far out as 1.5 virial radii for 
halos that are on first infall ( |Hahn et al.|2009| l. We find that first- 
infall centrals within 1.5 i?200m have average halo masses that are 
~ 10% lower than isolated centrals, so this does contribute to 
Fig. [9] but is not the dominant effect. We also checked that the 
distribution of stellar masses within our chosen range does not 
vary with distance to create additional signal. 

7 The choice of velocity cut represents a trade-off between pu- 
rity (intrinsic signal) and completeness (statistics), and the choice 
should be optimized for each observational measurement. 
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los with M2oom = 10 Mq, the number in each such dis- 
tance bin increases to ~ 5000. Fig. [9] (bottom) shows average 
halo masses for this much larger sample. While the average 
is 10 - 15 percentage points higher, the order-of-magnitude 
increase in galaxy count makes such a lensing measurement 
potentially feasible with existing SDSS data, and upcom- 
ing measurements from deeper imaging surveys such as the 
Dark Energy Survey (DES), Hyper Suprime-Cam (HSC), 
and the Large Synoptic Survey Telescope (LSST) should 
provide strong constraints on these halo masses of ejected 
satellites (e.g., Li et al.|2013 1. 

Indeed, galaxy-galaxy lensing measurements already 
show reduced subhalo masses for satellite galaxies within 



individual clusters (Limousin et al. 2007 Natarajan et al 



20091, and Gillis et al. (20131 measured reduced (sub)halo 



masses for galaxies in high-density environments (defined 
by local galaxy density) from the Canada-France-Hawaii 
Telescope Lensing Survey (CFHTLenS). In addition to mea- 
surements of subhalo masses within groups/cluster, such as 
these, we emphasize the utility of measuring galaxies near 
groups/clusters: while the (sub)halo masses of ejected satel- 
lites are not as tidally truncated as those of satellites at the 
core of their group/cluster, galaxies near groups/clusters can 
provide a cleaner lensing signal because their contrast with 
the background density field is much stronger. Overall, these 
results suggest promise in measuring (sub)halo stripping in 
ejected satellites via galaxy-galaxy lensing, and we will pur- 
sue further investigation in future work. 

In summary, ejected satellites have systematically bi- 
ased halo masses, typically only 40% as high as all cen- 
tral galaxies of the same stellar mass, leading to a S3 2.5 x 
higher M sta r/Mh a io ratio. This truncation strongly manifests 
itself in the average halo masses of central galaxies out to 
S3 2.5i?200m beyond groups/clusters and should be measur- 
able via galaxy-galaxy lensing, particularly with measure- 
ments from surveys such as DES, HSC, and LSST. 



7 SUMMARY AND DISCUSSION 
7.1 Summary of Results 

Using a galaxy group/cluster catalog from SDSS to decom- 
pose galaxies into 'centrals' and 'satellites', we examined the 
quiescent fraction for galaxies with M sta r = 10 9,7-10 ' 5 Mq 
that reside near groups/clusters with M2oom = 10 12 ~ 15 Mq. 
We postulated that the observed enhancement of the 
quiescent fraction in central galaxies beyond the virial 
radius of host halos is caused by ejected satellites, that is, 
galaxies that fell into a host halo and then orbited back out 
beyond the virial radius, becoming again central galaxies in 
their own halo. We used a cosmological iV-body simulation 
to examine their orbital histories and mass evolution, and 
we examined scenarios for their SFR evolution, which we 
tested by applying the same group-finding algorithm to the 
simulation to make robust comparisons with SDSS. Our 
primary results are as follows. 

Central galaxies near more massive host halos are more 
likely to be quiescent. Considering all galaxies regardless of 
type, an enhanced quiescent fraction (beyond the cosmic 
average at a given stellar mass) extends out to ~ 6i?200m 



around massive host halos. However, this enhancement is 
caused partially by satellites that reside in neighboring host 
halos and the fact that host halos are spatially clustered. 
Considering only central galaxies, which reside in their own 
host halo, an enhanced quiescent fraction extends out to 
S3 5i?200m and is especially strong within 2.5i?200m- 

A large fraction of central galaxies near massive host 
halos are ejected satellites. The orbits of ejected satellites ex- 
tend to S3 2.5i?200m beyond their host halo, though because 
of satellites ejected from neighboring host halos, ejected 
satellites represent a significant fraction of all central galax- 
ies out to much larger distances. Satellites of lower mass and 
those in/near more massive host halos are more likely to ex- 
perience an ejected phase, a likely result of the decreased 
efficiency of dynamical friction braking in those regimes. In 
the extreme case, ejected satellites represent 40% of all cen- 
tral galaxies of M s t ar < 10 10 Mq within 2.5 i?200m of clusters. 

The orbital histories of ejected satellites are consistent 
with satellite-like environmental quenching. After first in- 
fall, ejected satellites almost always (~ 90%) experienced a 
single pericentric passage, typically to ~ l/3i?200m, spend- 
ing roughly a virial crossing time, ~ 3 Gyr, within their 
host halo. Almost all remain on bound orbits and (will) re- 
enter the host halo ~ 3.5 Gyr after ejection. Their typical 
time since first infall is ~ 5 Gyr, longer than the satellite 
quenching timescale, and almost all (90%) continue to lose 
(sub)halo mass after ejection and do not re-accrete mass. 

The quenching of star formation in ejected satellites is 
consistent with occuring in the same manner as in normal 
satellites, and this can explain the enhanced quiescent frac- 
tion for central galaxies near massive host halos. We applied 
the same model for satellite quenching, based on time since 
first infall, that we developed for normal satellites in |Pa-| 
|per III| to ejected satellites. This simple extension of our 
model naturally and accurately accounts for the enhanced 
quiescent fraction in central galaxies out to S3 5i?200m- 

Effectively all environmental quenching is consistent 
with occurring within a host halo's virial radius. The suc- 
cess of our simple and natural extension of satellite quench- 
ing timescales to ejected satellites leaves little-to-no room 
for additional environmental effects for quenching central 
galaxies prior to first virial infall. 

Ejected satellites have significantly biased halo masses, 
and this is potentially observable. The (sub)halos of ejected 
satellites lose significant mass from tidal stripping after both 
first infall and ejection. To the extent that the stellar mass 
remains unstripped, ejected satellite have a S3 2.5 times 
higher M star /Mh a i ratio (or 60% lower halo mass) than all 
central galaxies of the same stellar mass. A reduction of 
average halo mass extends to central galaxies S3 2.5i?200m 
beyond massive host halos and is potentially observable via 
galaxy-galaxy weak lensing. 



7.2 Discussion 

7.2.1 Comparison with previous work 

Our results broadly agree with those of previous work on 
the SFRs of ejected satellites, though our approach repre- 
sents a significant advancement: using a larger galaxy sam- 
ple and examining a broader range of galaxy and host halo 
masses, accounting for correlated structure and neighboring 
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host halos through our SDSS group catalog, tracking indi- 
vidual satellite orbital histories using a simulation, apply- 
ing an empirically constrained model for SFR evolution and 
quenching in satellites, and robustly comparing our model 
results with SDSS through a simulation mock group catalog. 



Mahajan et al. (20111 pursued a related analysis of 



galaxies near massive clusters using the SDSS group cat- 
alog of Yang et al. (20071. By statistically comparing the 



phase-space distribution (in redshift space) of galaxies near 
clusters with that of particles taken from a smooth parti- 
cle hydrodynamics (SPH) simulation, they also argued that 
ejected satellites are quenched significantly relative to in- 
falling galaxies, and that star formation quenching occurs 
after a single pericentric passage. 

Y. | Wang et al.| ( |2009| > used the SDSS group catalog of 
Yang et al. ( 2007 1 to examine the g — r colors of even fainter 
galaxies (M r — 51og(/i) > —17). They found that the reddest 
population of such central galaxies lies within ~ 3i?i80m 
of a nearby massive host halo and therefore argued that 
they are likely ejected satellites. Curiously, they noted that 
more than half of their red, central galaxies do not lie within 
~ 3 -Ri80m of a more massive halo. However, the more recent 



analysis of Geha et al. (2012 I, using spectroscopic measures 



of star formation, found no instances of an isolated galaxy 
(with no more massive galaxy within 1.5 Mpc) at M sta r < 
10 9 Mq being quiescent. This suggests that the population 
of isolated but red galaxies in Y. Wang et al. is likely a result 
of dust reddening — a significant effect for low-mass galaxies 
(e.g., Mailer et al. 2009 1 — and not physical quenching. Thus, 
the results of |Geha et al.| suggest that any central galaxy at 
A/star < 10 9 Mq that is quiescent is an ejected satellite. 



Relatedly, |Teyssier et al. (20121 argued that satellite 
ejection likely explains the presence of several dwarf galax- 
ies that are beyond the virial radius of the Milky Way and 
have depleted HI masses and old stellar populations, based 
on these galaxies' observed phase-space coordinates. Our re- 
sults, for a much larger sample of more massive galaxies, 
firmly support this scenario. 



Knebe et al. (20111 examined the mass-to-light ratios 



of ejected (backsplash) versus infalling central galaxies in 
an SPH simulation of a Local Group analogue. They also 
found that ejected satellites have a higher ratio of stellar to 
halo mass than infalling central galaxies, but interestingly, 
they measured the mass-to-light ratio within the extent of 
the galactic stellar component (farthest star particle). Their 
result suggests that, in addition to stripping halo mass, tidal 
forces might heat the orbits of stars in the discs of ejected 
satellites, which could provide another means to differenti- 
ate observationally ejected satellites from infalling centrals 
based on stellar concentration and/or kinematics. 



Bahe et al. (2013) examined extended environmental 



processes on galaxies in the GIMIC suite of SPH simulations. 
While they found some stripping of extended gas (beyond 
the disk) in the halos of infalling galaxies prior to virial 
crossing from ram-pressure in filaments, the effect on star 
formation before virial crossing was negligible for galaxies in 
our mass range, in agreement with our results. They did see 
stronger reduction of star formation prior to virial crossing 
for galaxies of significantly lower mass, which could be tested 
by observations of lower mass galaxies. 



7.2.2 Robustness of results 

We have argued that ejected satellites, which passed within 
the virial radius of a more massive host halo, can account 
for essentially all galaxy quenching beyond the virial radius. 
However, we cannot entirely rule out environmental pro- 
cesses that extend beyond the virial radius, quenching cen- 
tral galaxies prior to crossing inside. Our model parametrizes 
the quenching of star formation in ejected satellites in terms 
of time since first infall (first virial radius crossing), with an 
initial 'delay' time over which star formation is unaffected, 
followed by rapid quenching. From |Paper III| this model suc- 
cessfully accounts for the SFR distribution of satellite galax- 
ies in groups/clusters as a function of both galaxy and host 
halo mass, and, as |Paper V| will show, successfully describes 
the dependence of the quiescent fraction on halo-centric dis- 
tance within groups/clusters. However, it may be possible 
that other models, in which the quenching process starts at 
a farther distance, work as well. We will test and compare 
other physical models for satellite quenching in Pap er V] 

None-the-less, our model, in which the environmental 
quenching process starts at first virial infall, can explain 
the SFRs of central galaxies beyond groups/clusters, and 
we argue that this is the simplest and most natural expla- 
nation for the extended environmental quenching, particu- 
larly given what we already know about satellite quench- 
ing within the virial radius. The SPH results of Bahe ct al. 
( |2013[ ) also support this picture for SFR in galaxies in our 
mass range. Thus, to the extent that the extension of our 
model to ejected satellites is correct, the agreement of our 
model with SDSS means that any environmental effects be- 
yond the virial radius are sub-dominant. 

Relatedly, we note that the exact choice of virial ra- 
dius is somewhat arbitrary, and the complications of ejected 
satellites could be 'defined' away, to some extent, by using 
a more liberal virial definition to encompass all associated 
satellites. However, this alternative has several drawbacks. 
Our virial definition of 200p m already is broader than many 
other commonly used ones, such as 200p c . Naively extrapo- 
lating an NFW density profile, a virial definition whose ra- 
dius extends to 2.5 -ffeoom — the extent of ejected satellites — 
would contain an average density of only ~ 20/9 m , a density 
more closely associated with filaments. A larger virial ra- 
dius also would encompass more galaxies that are not (yet) 
physically associated with the host halo and thus have not 
(yet) likely been affected environmentally. Thus, the choice 
of virial radius necessarily represents a trade-off between pu- 
rity and completeness of the associated satellite population, 
and a proper treatment of ejected satellites will be impor- 
tant for any reasonable virial definition. 

Additionally, while our group-finding algorithm as- 
sumes that host halos are spherical, their density contours 
and virial shock fronts can be ellipsoidal, meaning that the 
'physical' extent of a host halo can extend to more (or less) 
than a spherical virial radius. We leave investigation of al- 
ternative group finders, in the context of ejected satellites, 
to future work. 



7.2.3 Broader implications 

Our results highlight the importance of environmental his- 
tory in understanding not only star formation but also stel- 
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lar mass growth in galaxy evolution. Central galaxies that 
are really ejected satellites represent a challenge for the stan- 
dard 'halo model' and halo occupation distribution (HOD) 
formalism (see Cooray fc Sheth|2002 for review). Both the 
SFR/color and stellar mass of ejected satellites do not de- 
pend simply on their halo mass, but also on their forma- 
tion history. As Fig. [8] showed, ejected satellites have sig- 
nificantly higher stellar mass than their instantaneous halo 
mass would imply, and they are more likely to be quies- 
cent/red than normal central galaxies of the same stellar 
mass. Thus, any modeling approach that seeks to assign 
stellar mass or SFR/color to halos based only on their in- 
stantaneous virial mass, as in a standard halo-model/HOD 
approach, will underestimate both the stellar mass and the 
quiescent/red fraction for such central galaxies. Moreover, 
because less massive central galaxies are intrinsically less 
likely to be quiescent /red, both effects bias in the same di- 
rection. This bias will be most important on quasi-linear 
scales in the vicinity of massive host halos ('1-halo' to '2- 
halo' transition), where the ejected fraction is 10 - 40%. 
It may be possible to modify the standard halo model to 
incorporate analytically an associated satellite population 
that extends out to 2.5 virial radii and spatially overlaps 
(in a statistical sense) with the population of distinct host 
halos, though such a modification would require a careful 
recalibration of the mass function and bias of host halos. 
Overall, ejected satellites represent a population for which 
it is particularly important to follow the formation history 
of a (sub)halo to understand the galaxy inside. 

Relatedly, the biased M star — Mhaio relation for ejected 
satellites, if unaccounted for, could manifest itself as in- 
creased scatter in the M sta r — Mhaio relation if averaging 
over all central galaxies in halo-model/HOD-based analy- 
ses. While this effect is unlikely to be significant at M star > 
10 10 M.0, where the fraction of all central galaxies that are 
ejected satellites is < 5% (Fig.[3|, we expect it to be increas- 
ingly significant at lower Af sta r where the overall ejected 
fraction is closer to 10%. Note that if 5 - 10% of galaxies 
at a given stellar mass have systematically 60% lower halo 
mass, this alone translates to a scatter of 0.09 - 0.12 dex 
in averaging the population. This represents a non-trivial 
fraction of empirically measured values of 0.15 - 0.2 dex (see 



[ 2.2.4 1. Thus, while we do not expect ejected satellites to be 
the dominant source of scatter in the Af sta r — Afhaio relation, 
they may be a significant component if unaccounted for. 

Finally, we note the potential connection of these re- 
sults to 'galactic conformity', that the properties of central 
galaxies, such as SFR/color, are correlated with those of 
their satellites ( |Weinmann et al."| [2006), or more generally, 
with all neighboring galaxies out to ~ 4 Mpc (Kauffmann 
et al. 20131. As iKauffmann et al.l noted 



20131. As |Kauffmann et al.| noted, this conformity 
signal is most pronounced for less massive galaxies that 
have low SFRs, and can be stronger at larger separations 
(> 1 Mpc). Ejected satellites provide a natural population to 
help explain this effect, both because (1) they are relatively 
low-mass, quiescent central galaxies that reside near more 
massive halos, which host many quiescent satellite galaxies, 
and (2) if ejected satellite halos contain their own satellites, 
then both the central and satellite galaxies in such a halo 
would have passed through a more massive host halo, thus 
both galaxy types are more likely to have environmentally 



quenched (correlated) star formation. However, Kauffmann 



|et al.| noted that a full explanation of their observed con- 
formity effect would require more than half of all low-mass 
central galaxies to be ejected satellites, and our results sug- 
gest that this fraction is always < 10% for galaxies with 
A/ Star > 10 9 Mq . None-the-less, we fully expect that ejected 
satellites are a population that should exhibit some confor- 
mity signal and thus would provide some contribution to 
the overall conformity effect. We will investigate this issue 
in future work. 



ACKNOWLEDGMENTS 

We thank Michael Blanton, David Hogg, and collaborators 
for publicly releasing the NYU VAGC, Jarle Brinchmann 
and the MPA-JHU collaboration for publicly releasing their 
spectral reductions, and Martin White for simulation data. 
We thank Maria Geha and Kathryn Johnston for enlighten- 
ing conversations, and Maria Geha for comments on an early 
draft. The simulation was analyzed at the National Energy 
Research Scientific Computing Center. 



REFERENCES 

Abazajian K. N., et al., 2009, ApJS, 182, 543 
Bahe Y. M., McCarthy I. G, Balogh M. L., Font A. S., 
2013, MNRAS 

Balogh M. L., Morris S. L., Yee H. K. O, Carlberg R. G, 

Ellingson E., 1997, ApJ, 488, L75+ 
Balogh M. L., Navarro J. F., Morris S. L., 2000, ApJ, 540, 

113 

Berrier J. G, Bullock J. S., Barton E. J., Guenther H. D., 

Zentner A. R., Wechsler R. H., 2006, ApJ, 652, 56 
Blanton M. R., Berlind A. A., 2007, ApJ, 664, 791 
Blanton M. R., Moustakas J., 2009, ARA&A, 47, 159 
Blanton M. R., Roweis S., 2007, AJ, 133, 734 
Blanton M. R., Schlegel D. J., Strauss M. A., Brinkmann 
J., Finkbeiner D., Fukugita M., Gunn J. E., Hogg D. W., 
Ivezic Z., Knapp G. R., Lupton R. H., Munn J. A., Schnei- 
der D. P., Tegmark M., Zehavi I., 2005, AJ, 129, 2562 
Boylan-Kolchin M., Ma C.-R, Quataert E., 2008, MNRAS, 
383, 93 

Brinchmann J., Chariot S., White S. D. M., Tremonti G, 
Kauffmann G, Heckman T., Brinkmann J., 2004, MN- 
RAS, 351, 1151 

Cacciato M., van den Bosch F. G, More S., Mo H., Yang 
X., 2013, MNRAS, 430, 767 

Cameron E., 2011, PASA, 28, 128 

Chabrier G, 2003, PASP, 115, 763 

Conroy C, Wechsler R. H., Kravtsov A. V., 2006, ApJ, 
647, 201 

Cooray A., Sheth R., 2002, Phys. Rep., 372, 1 

Davis M., Efstathiou G, Frenk C. S., White S. D. M., 1985, 

ApJ, 292, 371 
Dressier A., 1980, ApJ, 236, 351 
Dressier A., Gunn J. E., 1983, ApJ, 270, 7 
Drory N., Bundy K., Leauthaud A., Scoville N., Capak P., 

Ilbert O., Kartaltepe J. S., Kneib J. P., McCracken H. J., 

Salvato M., Sanders D. B., Thompson D., Willott C. J., 

2009, ApJ, 707, 1595 



© 2013 RAS, MNRAS 000,[T|fl4l 



14 Wetzel, Tinker, Conroy & van den Bosch 



Fraternali F., Tolstoy E., Irwin M. J., Cole A. A., 2009, 

A&A, 499, 121 
Geha M., Blanton M. R., Yan R., Tinker J. L., 2012, ApJ, 

757, 85 

Gill S. P. D., Knebe A., Gibson B. K., 2005, MNRAS, 356, 
1327 

Gillis B. R., Hudson M. J., Erben T., et al., 2013, MNRAS 

Hahn O., Porciani C, Dekel A., Carollo C. M., 2009, MN- 
RAS, 398, 1742 

Hansen S. M., Sheldon E. S., Wechsler R. H., Koester B. P., 
2009, ApJ, 699, 1333 

Hogg D. W., Blanton M. R., Brinchmann J., Eisenstein 
D. J., Schlegel D. J., Gunn J. E., McKay T. A., Rix H.- 
W., Bahcall N. A., Brinkmann J., Meiksin A., 2004, ApJ, 
601, L29 

Jiang C. Y., Jing Y. P., Faltenbacher A., Lin W. P., Li C, 

2008, ApJ, 675, 1095 

Kauffmann G., Li C, Zhang W., Weinmann S., 2013, MN- 
RAS, 430, 1447 
Kauffmann G., White S. D. M., Heckman T. M., Menard 

B. , Brinchmann J., Chariot S., Tremonti C, Brinkmann 
J., 2004, MNRAS, 353, 713 

Knebe A., Libeskind N. I., Knollmann S. R., Martinez- 
Vaquero L. A., Yepes G., Gottlober S., Hoffman Y., 2011, 
MNRAS, 412, 529 

Leauthaud A., Tinker J., Bundy K., et al., 2012, ApJ, 744, 
159 

Li C, White S. D. M., 2009, MNRAS, 398, 2177 
Li R., Mo H. J., Fan Z., Yang X., van den Bosch F. C, 
2013, MNRAS 

Limousin M., Kneib J. P., Bardeau S., Natarajan P., Czoske 
O., Smail I., Ebeling H., Smith G. P., 2007, A&A, 461, 881 
Ludlow A. D., Navarro J. F., Springel V., Jenkins A., Frenk 

C. S., Helmi A., 2009, ApJ, 692, 931 

Mahajan S., Mamon G. A., Raychaudhury S., 2011, MN- 
RAS, 416, 2882 
Mailer A. H., Berlind A. A., Blanton M. R., Hogg D. W., 

2009, ApJ, 691, 394 

Mamon G. A., Sanchis T., Salvador-Sole E., Solanes J. M., 

2004, A&A, 414, 445 
Mandelbaum R., Seljak U., Kauffmann G., Hirata C. M., 

Brinkmann J., 2006, MNRAS, 368, 715 
Mateo M. L., 1998, ARA&A, 36, 435 

More S., van den Bosch F. C, Cacciato M., Mo H. J., Yang 

X., Li R., 2009, MNRAS, 392, 801 
Natarajan P., Kneib J., Smail I., Treu T., Ellis R., Moran 

S., Limousin M., Czoske O., 2009, ApJ, 693, 970 
Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ, 490, 

493 

Noeske K. G., Weiner B. J., Faber et al., 2007, ApJ, 660, 
L43 

Oemler Jr. A., 1974, ApJ, 194, 1 

Reddick R. M., Wechsler R. H., Tinker J. L., Behroozi P. S., 

2012, ArXiv e-prints 
Sales L. V., Navarro J. F., Abadi M. G., Steinmetz M., 

2007, MNRAS, 379, 1475 
Solanes J. M., Sanchis T., Salvador-Sole E., Giovanelli R., 

Haynes M. P., 2002, AJ, 124, 2440 
Teyssier M., Johnston K. V., Kuhlen M., 2012, MNRAS, 

426, 1808 

Tinker J., Kravtsov A. V., Klypin A., Abazajian K., War- 
ren M., Yepes G., Gottlober S., Holz D. E., 2008, ApJ, 



688, 709 

Tinker J., Wetzel A., Conroy C, 2011, ArXiv e-prints 
Trujillo-Gomez S., Klypin A., Primack J., Romanowsky 

A. J., 2011, ApJ, 742, 16 
Vale A., Ostriker J. P., 2006, MNRAS, 371, 1173 
van den Bosch F. C, Aquino D., Yang X., Mo H. J., 

Pasquali A., Mcintosh D. H., Weinmann S. M., Kang X., 

2008, MNRAS, 387, 79 
von der Linden A., Wild V., Kauffmann G., White S. D. M., 

Weinmann S., 2010, MNRAS, 404, 1231 
Wang H., Mo H. J., Jing Y. P., 2009, MNRAS, 396, 2249 
Wang Y., Yang X., Mo H. J., van den Bosch F. C., Katz 

N., Pasquali A., Mcintosh D. H., Weinmann S. M., 2009, 

ApJ, 697, 247 

Warnick K., Knebe A., Power C, 2008, MNRAS, 385, 1859 
Weinmann S. M., van den Bosch F. C, Yang X., Mo H. J., 

2006, MNRAS, 366, 2 
Wetzel A. R., 2011, MNRAS, 412, 49 

Wetzel A. R., Tinker J. L., Conroy C, 2012, MNRAS, 424, 
232 

Wetzel A. R., Tinker J. L., Conroy C, van den Bosch F. C, 

2012, MNRAS, p. submitted 

Wetzel A. R., Tinker J. L., Conroy C, van den Bosch F. C, 

2013, MNRAS 

Wetzel A. R., White M., 2010, MNRAS, 403, 1072 
White M., 2002, ApJS, 143, 241 

White M., Cohn J. D., Smit R., 2010, MNRAS, 408, 1818 
Yang X., Mo H. J., van den Bosch F. C, 2008, ApJ, 676, 
248 

Yang X., Mo H. J., van den Bosch F. C, 2009, ApJ, 693, 
830 

Yang X., Mo H. J., van den Bosch F. C, Jing Y. P., 2005, 

MNRAS, 356, 1293 
Yang X., Mo H. J., van den Bosch F. C, Pasquali A., Li 

C, Barden M., 2007, ApJ, 671, 153 



© 2013 RAS, MNRAS 000,[T]{T4] 



